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Abstract 

A generic formalism is presented for the time-dependent or time-integrated decays of 
any coherent p°p° system (P° = K° , D° , B®, or B®). To meet various possible measure- 
ments at asymmetric B factories, we reanalyze some typical signals of CP violation in 
the coherent B®B® transitions. The advantage of proper time cuts is illustrated for mea- 
suring mixing parameters and CP violation. We show that the direct and indirect CP 
asymmetries are distinguishable in neutral B decays to CP eigenstates. The possibility 
to detect the CP- forbidden processes at the T(45) resonance is explored in some detail. 
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1. Introduction 



It is known in particle physics that mixing between a neutral meson P° and its CP-conjugate 
state P° provides a mechanism whereby interference in the decay amplitudes can occur, leading 
to the possibility of CP violation [1]. To date, K° — K° mixing and B d — B d mixing have been 
measured, and the CP violating signal induced by K° — K° mixing has been unambiguously 
established [2]. Compared with B d — B d mixing, B° s — B° s (D° — D°) mixing is expected to be 
quite large (very small) in the context of the standard electroweak model [3]. Today the B d — B d 
system is playing an important role in studying flavour mixing and CP violation beyond the 
neutral kaon system. The B° s — B® and D° — D° systems are more interesting in practice for 
probing new physics that is out of reach of the standard model predictions. 

A feasible way to study CP violation arising from P° — P° mixing is to measure the coherent 
decays of P°P° pairs produced at appropriate resonances. For instance, 

(f) -> K°K° , V" -> D°D° , T(4S) -> B° d B° d , T(5S) -> B° S B° S . 

At present, efforts are underway to develop asymmetric B factories at KEK and SLAC lab- 
oratories [4], while asymmetric <fi factory options are also under consideration [5]. The main 
purpose of these machines is to probe CP violation (and to test other discrete symmetries or 
conservation laws) by measuring the time-dependent transitions. By now some phenomenolog- 
ical analyses of coherent K°K° and B d B d decays have been made in the literature [6-8]. These 
works have outlined the main features of CP violation in the K° — K° or B d — B d systems, 
although many of their formulae and results rely on the characteristics of the system itself or 
some model-dependent approximations. A generic formalism, which can describe the common 
properties of coherent P°P° decays, is still lacking. In addition, little attention has been paid 
to the advantage of proper time cuts for measuring the mixing parameters and CP asymmetries 
in coherent weak decays. 

In this paper we shall present a generic and concise formalism for the time-dependent or 
time-integrated decays of any coherent P° — P° system. This formalism should be very useful 
for phenomenological applications, because it is independent of the contexts of the standard 
model or its various non-standard extensions. To meet various possible measurements at the 
forthcoming B factories, we shall carry out a reanalysis of the typical signals of CP violation 
manifesting in B d — B d mixing, in B d vs B d decays to CP eigenstates, and in CP- forbidden 
transitions at the T(4S) resonance. Our analysis differs from the previous ones in the following 
three aspects: (1) we illustrate the advantage of proper time cuts for measurements of the 
mixing parameters and CP violating asymmetries; (2) we highlight the distinguishable effect 
of direct CP violation on CP asymmetries in neutral B decays; and (3) we explore in some 
detail the possibility to detect the CP-forbidden decays at B factories. This work concentrates 
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on analytical studies. A more comprehensive discussion about the present topic, together with 
numerical (model-dependent) predictions, will be given elsewhere. 



2. Formalism of coherent P°P° decays 



(i) 



The time-dependent wave function for a Pphy S Pphys P a i r a ^ res ^ can be written as 

[|P p ° hys (K,0> ® |P p ° hys (-K,t)) + C|P p ° hys (-K,*)> ® |P p ° hys (K,t)); , 

where K is the three-momentum vector of the P meson, and C = — or + is the charge- 
conjugation parity of the Pph ys Pphys pair. The proper time evolution of an initially (t = 0) pure 
P° or P° is given by 



^hysW) = g + (t)\P°) + (q/p)g-(t)\P }, 



IPpVsW) = (p/ q )g-(t)\P°) + g+(t)\P ), 



(2) 



(=)„ 



where the mixing parameters p and q connect the flavour eigenstates | P°) to the mass eigen- 
states |Pi, 2 ) through |Pi) = p\P°) + q\P°) and |P 2 ) = p\P°) - q\P°); and 



g ± ( t ) = -e-i^y e +( tAm -^)I ±e -(iAm-^) 
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(3) 



Here we have defined m = (mi+m2)/2, Am = (m2-mi), T = (ri + r2)/2, and Ar = (ri — r 2 ), 
where I\2 and m^2 are the widths and masses of Pi )2 ■ 

Now we consider the case that one of the two P mesons (with the momentum K) decays to 
a final state fi at proper time t\ and the other (with — K) to f 2 at t 2 . j\ and f 2 may be either 
hadronic or semileptonic states. After a lengthy calculation, the joint decay rate for having 
such an event is given as 



R(/i,*i;/2,f 2 )c oc \A fl \ 2 \A h \ 2 e- n + 

Kcl 2 - ICd 2 ) cos(Amt c ) + 21m (C c Cc) sin(Amt c ) 



where 



and 



a u = (miP ) , a /4 = (m\p°) , p h = -r > (< = 1.2) 



(4) 

(5) 
(6) 



*c = t 2 + Cfi, £c = (p/q) + C(q/p)p fl p h (c = Pf 2 + C Ph- 

The time-independent decay rate is obtainable from Eq. (4) by integrating R(/i, ti, f 2 , t 2 )c 
over ti and t 2 : 



R(h,f 2 )c oc I^PI^I 2 
1-Cx 2 



\Zc + (c\ 2 \Zc-(c\ 2 



+ 



(1 + x 2 ) 2 



2(l + y)(l + Cy) 2(l-y)(l-Cy) 
(l^| 2 -|Cc| 2 ) + ^^Im(^Cc) 



(7) 
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where x = Am/T and y = Ar/(2r) are two measurables of the P° — P° system. Note that 
Eqs. (4) and (7) are useful at both symmetric and asymmetric flavour factories. 

An asymmetric e + e~ collider running at the threshold of production of (-Pph ys -Pph y s)c pairs 
will offer the possibility to measure the decay-time difference t- = (t 2 — h) between Pp hys — * fi 
and -Pphys — >• /2 - It is usually difficult to measure the t + = (t 2 + ti) distribution in either linacs 
or storage rings, unless the bunch lengths are much shorter than the decay lengths [4,5,9]. Hence 
it is more practical to study the t- distribution of the joint decay rates. Here and hereafter 
we use t to denote t- for simplicity. Integrating R(/i,ti; /2,t 2 )c over t+, we obtain the decay 
rates (for C — ±) as follows: 



R(/i,/ 2 ;i)- oc 



l^/i! 2 |^/ 2 | 2 e 



r|t| 

IC-I 2 



C- 



\\H- + t-\ 2 e- yVt + \\li 
cos(xH) + 21m (C(-) sin(a;rt) 



(8) 



and 



R(/i,/ 2 ;*)+ oc 



A /2 | 2 e- r l*l 



ie+i 2 -ic + i 2 



2(1 + 2/) 



+ 1 p -»r|t| 



ie + -c + p 

2(1 -y) 



3 +yT\t\ 



21m U* + C, 

cos (xT\t\ + <f) x ) H , ' %m.{xY\t\ + (f> x ) 



(9) 



where <p x = arctana; signifies a phase shift. One can check that integrating K(fi, f 2 ',t)c over 
t, where t G (— oo,+oo), will lead to the time-independent decay rates R(fi,f 2 )c m Eq. (7). 
Eqs. (8) and (9) are two basic formulae for investigating coherent B°B° (or K°K°) decays at 
asymmetric B (or 0) factories. 

Another possibility is to measure the time-integrated decay rates of (PphysPphys)^ with a 
proper time cut, which can sometimes increase the sizes of CP asymmetries [10]. In practice, 
appropriate time cuts can also suppress background and improve statistic accuracy of signals. 
If the decay events in the time region t G [-No, +oo) or t G (— oo, —to] are used, where to > 0, 
the respective decay rates can be defined by 



R(/i,/ 2 ;+fo)c 
R(/i, / 2 ; —to)c 

From Eqs. (8) and (9) we obtain 



RC^/ajito)- oc \A fl \ 2 \A f2 \ 2 e- rt ° 



ie-i 2 -ic-i 



+t 



-to 



R(hJ 2 ;t)cdt, 
R(h,f 2 ;t) c dt. 



4(1 + y) 



; e -yTt 0+ \Z-T(- : ( 



2 V / TT 



cos (xFto + <f>x) ± 



4(1-2/) 

im (ec-) 



(10) 



x/ 



sin(a;rto + 4> x ) 



;n) 
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and 



R(/i,/ 2 ;±*o)+ « \A fl \ 2 \A 



fit 



2 e -rt 



ie + + c + i 2 

4(1 + y) 2 



ie + i 2 -ic + i 

2(1 + x 2 ) 



cos (xTt + 2(j) x ) + 



ie + -c + i 2 

' 4(1 -yf 



1 + x 2 



sin(xrt 



(12) 



It is easy to check that 



R(/i,/ 2 ;+o)c + R(/i,/ 2 ;-o) c = R(A,/ : 



2)C 



(13) 



One can observe that in R(/i, A> ±^o)c different terms are sensitive to the time cut to i n 
different ways. Thus it is possible to enhance a CP violating term (and suppress the others) 
via a suitable cut to ■ 

The formulae given above are applicable to all coherent decays of the K° — K°, D° — D°, 



B® — B®, and B® — B° s systems. Denoting the decay amplitudes of P n —> /j by Ai 



(») 



n, i 



1,2) 



and the ratio of to A K j~.' by rjf, : one can also express the joint decay rates in terms of Aj 
and rjj. through the following transformations: 



i(D 



(n) 



1 

2p 



Af + A$ 



Ai 



1 

2^ 



.a) _ A m 

1 1 f y l jc 

J i 



Pfi 



p 



(14) 



9 1 + ^ 

Such notations are usually favoured in the — A' system [6]. In the following we shall apply 
the above formalism to the coherent (B^ phys B^ phys )c decays and CP violation at the T(45) 
resonance, a basis of the forthcoming B factories [4,9]. 



3. Signals of CP violation at B factories 

The unique experimental advantages of studying 6-quark physics at the T(4S) resonance 
are well known. For symmetric e + e~ collisions the produced B mesons are almost at rest and 
their mean decay length is only about 20/xm, a distance which is insufficient for identifying the 
decay vertices or measuring the decay time difference [9,10]. If the colliding e + and e~ beams 
have different energies, the product of collisions will move with a significant relativistic boost 
factor in the laboratory^ (along the direction of the more energetic beam). This can cause the 
two B mesons far apart in space, such that the distance between their decay vertices becomes 
measurable. It is then possible to study the time distribution of the joint decay rates and CP 
asymmetries. 

A. CP violation in B® — mixing 

^or a moving T(4S') system, the momentum of the B mesons in the T(4«S) rest frame can be ignored. This 
safe approximation has been discussed in Refs. [9,10]. 
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We first consider the joint decays (B dphys B dphys )c — > (I X^)(l Xjf), which lead to dilepton 
events in the final states. Keeping the AQ = AB rule and CPT symmetry, we have 



(l+Xf\H\B° d ) = (l-X+\H\B° d ) = A u ■ 
(l-X+\H\B°) = (l+Xr\H\B°) = 0, 



(15) 



where i = a or b. Subsequently we use N^ ± (t) and N£ - (i) to denote the time-dependent like- 
sign and opposite-sign dilepton numbers, respectvely. Similarly, let (±t ) and N^~(±t ) 
denote the time-integrated dilepton events with the time cut to • With the help of Eqs. (8) 
and (9), we obtain 



Nl+(t) ex \ P /q\ 2 \A la \ 2 \A lb \ 2 e- m [cosh(yTt)-cos(xTt)] 
N—(t) oc \q/p\ 2 \A la \ 2 \A lb \ 2 e- m [cosh(yrt) - cos(a?rt)] 
Ni~(t) oc 2\A la \ 2 \A lb \ 2 e- m [cosh{yTt) + cos{xTt)] ; 



(16) 



and 



N+ + (t) oc \p/q\ 2 \A la \ 2 \A lb \ 2 e-^ 
N-~(t) oc \q/p\ 2 \A la \ 2 \A lb \ 2 e^ 
N+-(t) oc 2\A la \ 2 \A lb \ 2 e- T ^ 



cosh(yT\t\) + y smh(yT\t\) cos(xT\t\ 



l-y 2 VTTxZ 
cos%r|t|) +ysinh(yr|t|) coa{xT\t\ + <f> x ) 



i-y 2 



VT+¥ 2 



(17) 



cosh(yT\t\) + y smh(yT\t\) + cos(xr|t| + <f> x ) 



l-y* 



If y is not very small in comparison with x, its size and sign should be (in principle) determinable 
from the above two equations!. 

Since both N^(t) and N£~(t) are even functions of proper time i, one finds 



N 



c y^o 



c^-to) = iHHto) , ^-(+to) = N+~Ho) = ~N+-(i ) • (18) 



The time-integrated observables of CP violation and B® — B® mixing can be defined as: 



c N+ + ((„)+Na-(i„) 

Using Eqs. (11) and (12) we obtain 

At- (to) = At"(*o) 



^" (to) = 



N+ + (t )+N c -(t ) 
Nj"(to) 



(19) 



b/g| 2 - lg/gf 

\p/q\ 2 + \q/p\ 



2 ' 



(20) 



2 A more detailed discussion has been given by Dass and Sarma in Ref. [7], where only the case of C 
was taken into account. 
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which is independent of the time cut to • The nonvanishing A c (to) implies CP violation in 
B® — B® mixing. In addition, we find 



SI' (to) = 



Ip/qI 2 + \q/p\ 2 cosh(yTt ) + ysmh(yTt ) - z cos(xTt + <p x ) 

2 cosh(yrt ) + y sinh(yn ) + z cos(xrt + 4> x ) ' foi \ 

\p/q\ 2 + \q/p\ 2 (1 + V 2 ) cosh(yrto) + 2ysmh( y rt ) - z 2 cos(xTt + 20.) 1 ] 

2 (1 + y 2 ) cosh(?/rto) + 2y sinh(yrt ) + z 2 cos(aTt + 2<p x ) ' 



where z — (1 — y 2 ) / Vl + x 2 . 

In the context of the standard model, \q/p\ 
Thus Eq. (21) is simplified as 



1 and y are two good approximations [3] . 



SI- (t ) 



Vl + x 2 - cos(xTt + (j) x ) to=0 

\/i + x 2 + cos(xrt + 4>x) 

(1 +x 2 ) - cos(xTt + 2(f) 



X 



t =0 



2 + x 2 ' 
3x 2 + x 4 

2 + x 2 + x 4 



(22) 



(1 + x 2 ) + cos(xTt + 2<t> x ) 

We show the evolution of S^ - (to) with to in Fig. 1, where the experimental input is x ~ 0.7 
[2]. One observes that an appropriate time cut can significantly increase the ratio of the same- 
sign dilepton events to the opposite-sign ones. Practically time cuts should be a useful way to 
enhance the signals of B% — B® mixing, only if the cost of the total number of events is not too 
large. 

B. CP asymmetries in B® vs B® decays to CP eigenstates 

Neutral B decays to CP eigenstates are favoured in both theory and experiments to study 
quark mixing and CP violation. At the T(4S) resonance, the produced B® phys and -B^ phys 
mesons exist in a coherent state untill one of them decays. Thus one can use the semileptonic 
decay of one B d meson to tag the flavour of the other meson decaying to a flavour-nonspecific 
hadron state. Let us consider the joint transitions (B® phys B® phys ) c (PX^f, where / 
denotes a hadronic CP eigenstate such as J/tpK s , D + D~ , or 7r + 7r~. To a good degree of 
accuracy in the standard model, we have \q/p\ ~ 1 and y xs 0. With the help of Eqs. (8) and 
(9), the time-dependent decay rates are given as 



R(^,/;0- oc \Ai 



\A f \ 2 e- T W 



1 + IA/l 2 . 1 



± 



cos(xTt) T- ImA/sin(a;rt) 



(23) 



and 



R(F,f;t) + « |^| 2 |A/| 2 e- r l«" 

1 



1 + l^fl 2 1 1 - |A f | 2 , , 

' " ± , COs(xT\t\ + <p x ) 



=F 



X/ 



2 "~ t/TTx^ 2 2 
ImA/ sin(a;r|t| + <p x ) 



(24) 
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where A/ = (q/p)pf ■ The time-dependent CP asymmetries, defined by 

R(l-J;t)c-R(l + ,f;t) c 



A c (t) 

can be explicitly expressed as 
A-(t) 



R(l~J;t)c + R(l + J;t) c ' 



(25) 



U f cos(xTt) + V/sin(a;rt) , 
A+(t) = 1 [U f cos(xT\t\ + + V/ sin(xr|t| + 



where 



1 ~ IA, 
1 + IA, 



-21mA/ 



(26) 



(27) 



1 + IA/l 2 

We find that Ac(t) contains both the direct CP asymmetry in the decay amplitude (|A/| ^ 1 
or Uf 7^ 0) and the indirect one from interference of mixing and decay (ImA/ ^ or V/ ^ 0). 
Measuring the time distribution of A±(t) can distinguish between these two sources of CP 
violation [8]. 

There are two ways to combine the time-integrated decay events (with the time cuts ±i ), 
leading to two types of CP asymmetries: 



A%\t ) 

With the help of Eqs. (11) and (12), we obtain 
(i) cos(xrto + 1 



R(l-J;+t ) + R(l-J; 


-to) 




R(l+,f;+to) + R(l + ,f; 


-h) 


R(l-,f;+t ) + R(l-,f; 


-to) 


+ 


R(l + ,f;+t ) + R(l + ,f; 


-to) 


R(l-,f;+t ) + R(l + ,f; 


-to) 




R(l+,f;+t ) + R(l-,f; 


-to) 


R(z- /;+<<,) + R(i+,/; 


-to) 


+ 


R(Z+,/;+* ) + R(Z-,/; 


-to) 



(28) 



A^'(t ) 



\fiTx 2 



-u 



A m {t Q ) 



sin(a:rto + (f) x ) 

VTTx 2 



V 



/ > 



and 



4p (to) = 



cos(xrt + 20 x ) sin(xno + 2<j) x ) 
1 + x 2 f 1 + x 2 



Vf, A? (t ) = 0. 



(29) 



(30) 



Clearly the asymmetries A- (to) and A- (to) signify direct and indirect CP violation, respec- 
tively. They can be separated from each other on the T(4S) resonance. In Fig. 2 we show 
the ratios A^\to) / dl\o) (n = 1 and 2) as functions of to • A proper time cut can certainly 
increase the CP asymmetries (at some cost of decay events). In practice, it can also suppress 
background and improve statistic accuracy of signals. Note that a suitable cut of decay time is 
(in principle) able to isolate the direct or indirect CP asymmetry in A+\to). For example, 

2 <Ar\ 1 7/ ,(1) ( 7T 2<p x \ 1 



7T 



xT xT 



1+x 2 



f » 



A 



2xT xF 



1 + x 2 



(31) 



At symmetric B factories, it is possible to measure a large CP asymmetry .4.^(0) 
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C. CP-forbidden transitions 



We finally consider the CP-forbidden decay modes 

(- B °,phys^d,phys)T ~^ (fafb)± , (32) 

where f a ,b denote the CP eigenstates with the same or opposite CP parities. It should be 
emphasized that for such decays CP violating signals can be established by measuring the 
joint decay rates other than the decay rate asymmetries [3]. In practice, this implies that 
neither flavour tagging nor time-dependent measurements are necessary. 

On the T(AS) resonance, the typical CP-forbidden channels include (P° iP hy S -^°,phys)- — * 
(//)+ with / = J/ipKsi J/i>K L) D + D~, and it + it~. Taking / = X Z cK s for example, where 
X 5c denote all the possible charmonium states that can form the odd CP eigenstates with 
K s (see Fig. 3), we have the safe approximation P x ^ c k s ~ — 1- m addition, we take y ~ 
and q/p e~ 2l/3 , where (3 corresponds to an inner angle of the Kobayashi-Maskawa unitarity 
triangle [2]. With the help of Eqs. (7) and (8), one obtains the branching fractions 

B(X 5c K s ,X 5c K s ;t)^ oc B 2 XecKs sm 2 (2(3) e - r ^ [1 - cos^H)] , 



B(X 5c K s ,X 5c K s )^ oc B^ cXs sin 2 (2^) T 



x 2 (33) 



x 2 



where Bx^ c k s denotes the branching ratio of B% — > X 5c Ks ■ Clearly the above joint decay rates 
are forbidden by CP symmetry ((3 = or ±7r). In practice, summing over the available final 
states X Sc K s can statistically increase the number of decay events: 

2 

B = Y,V(X-ccKs,X- cc K s )_ oc -J—sin 2 ^)^ 2 ^) . (34) 

cc ^cc 

Just above the T(4S f ) resonance, an interesting type of CP-forbidden channels should be 
(P° phys P°,phys)+ — *■ [{XccKs){XccKl)\-. Neglecting CP violation in the kaon system, we have 
Px- cc k l ~ ~Px- cc K s ~ 1 and B x, c Ks ~ b x Sc k l to a good degree of accuracy. From Eqs. (7) and 
(9), it is staight forward to obtain 

cos(xr|t| + 4> x 



B(X- cc K s ,X- cc K L ;t) + oc B^ s sin 2 (2/3)e- r l*l 1- 
B(X- CC K S ,X- CC K L ) + oc B M s sin 2 (2/3) T |^|^_ 
Summing over the possible states (X 5c Ks)(XccKl), we find 



\ZTTx 2 



(35) 



B + ee Y,V{X-ccK s ,X- cc K L ) + oc r ^^ I sin 2 (2/3)E(BL^) ■ ( 36 ) 



In Fig. 4 we show the relative sizes of the effective branching fractions B_ and B + in the region 
0.17 < sin(2/3) < 0.99, limited by the current data [11]. For our purpose, the suitable X gc states 
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include J/ip, ip' , ip" , r} c , r}' c: etcP. Since such channels occur through the same tree- level quark 
diagram (see Fig. 3), their branching ratios B x ^ c k s ( or ^x 5c k l ) are expected to be the same 
order. Thus a combination of the possible decays (B® phys B^ phys )c — > X 5c Ks or X 5c Kl should 
increase the decay events of a single mode by several times. In a similar way one can study 
the joint transitions (B^ phys B® phys ) c — > X ec K S 7r° or X 5 cK l ti . For a more detailed discussion 
about CP violation in the semi-inclusive decays B® vs B® — > (cc)Ks or (cc)Kl, we refer the 
reader to Ref. [12]. 

4. Summary 

In keeping with the experimental efforts to study flavour mixing and CP violation, we have 
presented a generic formalism for the time-dependent and time-integrated decays of all possible 
P° — P° systems. This formalism is useful for various phenomenological applications at the 
forthcoming flavour factories, where a large amount of coherent P°P° events will be produced 
and accumulated. In our calculations, the AQ = AP rule and CPT symmetry have been 
assumed. Relaxing these two constraints one can obtain the more general formulae, which 
should be useful for searching for CPT violation or AQ = AP breaking in the K° — K° [6] and 
B° - B° systems [13,14]. 

To meet various possible measurements of CP violation at asymmetric B factories, we have 
carried out a reanalysis of three types of CP violating signals manifesting in the coherent B^B^ 
transitions at the T(4S) resonance. Although some comprehensive works have been done on 
this topic, our present one differs from them in several aspects. We illustrate the advantage 
of proper time cuts for measuring the B^ — B® mixing parameters and CP asymmetries. It 
is shown that direct and indirect CP violating effects are distinguishable in neutral B decays 
to CP eigenstates. In addition, we explore the possibility to measure some CP-forbidden 
processes as a direct test of CP symmetry breaking at B factories. 

A more detailed study of the topic under discussion, together with some numerical predic- 
tions or estimates, is in preparation. 
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Figure 1: Ratios of the same-sign dilepton events to the opposite-sign ones as functions of the 
time cut t at the T(AS) resonance. 




Figure 2: Ratios of the CP asymmetries A- (to) to A- (0) as functions of the time cut t on 
the T(AS) resonance. 
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Figure 3: Quark diagrams for B® versus B® decays to X 5c Ks,l- Here X 5c denote the possible 
charmonium states that can form the odd (even) CP eigenstates with Ks (Kl)- 
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Figure 4: Relative sizes of the effective branching fractions Be (in arbitrary units) as functions 
of the CP violating angle (3 at the T(4S) resonance. 
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